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We analyze spectral properties of the ultrarelativistic (Cauchy) operator | A | 1 ^ 2 , provided its 
action is constrained exclusively to the interior of the interval [—1,1] C R. To this end both 
analytic and numerical methods are employed. New high-accuracy spectral data are obtained. A 
direct analytic proof is given that trigonometric functions cos(nnx/2) and sin(n7ra;), for integer n 
are not the eigenfunctions of |A|)( 2 , D = ( — 1,1). This clearly demonstrates that the traditional 
Fourier multiplier representation of lA^ 2 becomes defective, while passing from R to a bounded 
spatial domain D C R. 


I. INTRODUCTION 

Fractional (Levy-type) operators are known to be spatially nonlocal. This becomes an issue if confronted with a 
priori imposed exterior Dirichlet boundary data, which set the familiar (quantum) infinite well enclosure. Standard 
fractional Laplacians, at each instant of time, extend their nonlocal action to the entire real axis R and this property 
needs to be reconciled with the finite support D = (—1,1) C R of the infinite well with width equal 2. 

One of the obvious obstacles arising here is rooted in the fact, that the traditional Fourier multiplier representation 
of Levy operators is no longer operational in the finite interval, 0-0, see specifically HI- Compare e.g. also a 
discussion of that issue for the familiar (Laplacian-generated) quantum mechanical infinite well problem, see e.g. m 
and 0. 

To elucidate the above point, let us recall that the Fourier integral -A= J R Ik^ f {k)e~ lkx dk = ~d IM f{x)/d\x\^ = 

|A| M / 2 /(a;), where / stands for a Fourier transform of / £ L 2 (R) and g(k) = |fc| M /(fc) is presumed to be L 2 (R)~ 
integrable. It is |fc| M which plays the role of the pertinent Fourier multiplier. 

The above Fourier formula is quite often interpreted as a universal definition of both the fractional operator — | A| 
and that of the fractional derivative of the /r-th order, d fJ ,f(x)/d\x\ l ‘ = —|A|^/ 2 /(x), for ft £ (0,2). However, this 
definition is unquestionably valid only, if the fractional operator is defined on the whole real line R. More than that, 
it appears to be merely a specific admissible choice in the family of equivalent (while on R) definitions, 0- 

The nonlocal operator — lA^/ 2 , is known to generate two versions of so-called fractional dynamics (dimensional 
constants being scaled away): (i) the semigroup exp(— 1\ A|^/ 2 ) / and (ii) unitary exp(— it\ A|^/ 2 ) / ones. 

Apart from the unperturbed (free) case, the Fourier (multiplier) representation of the fractional dynamics has 
proved useful if an infinite or periodic support is admitted for functions in the domain, [s]. For the simplest quadratic 
(~ x 2 ) perturbation of the fractional Laplacian (the fractional oscillator problem), a complete analytic solution has 
been found in the ultrarelativistic (Cauchy) oscillator case, by resorting to Fourier space (and Fourier multiplier) 
methods. 

For more complicated perturbations, and likewise for a deceivingly simple problem of the fractional Laplacian in 
a bounded (spatial) domain, standard Fourier techniques have been found to be of a doubtful or limited use, 
Therefore, to keep spatial constraints under control, we turn over to a fully-fledged spatially nonlocal definition of the 
fractional Laplacian, that is well known in the mathematical and statistical physics literature 1-0 , while seldom 
invoked by quantum theory practitioners, see however [l] and references there in. 

Dating back to the classic papers 0,0, one interprets the fractional Laplacian — |A| M / 2 , g, £ (0,2) as a pseudo¬ 
differential (integral) spatially nonlocal operator and its action on a function from the L 2 (R) domain is commonly 
defined by employing the Cauchy principal value of the involved integral (evaluated relative to singular points of 
integrands) 


, 2/w = _ r(„+l) S in(,„/2) f /(;)- /M fc ( 1 ) 

7T J R 

For a rationale and a broader discussion of the uses (and misuses) of this formula, including its Fourier multiplier 
version, see 0. 

By departing from the general spatially nonlocal definition 0 we shall pass to the specialized Cauchy case (fi = 1) 
of the fractional Laplacian and next focus our attention on its properties under exterior Dirichlet boundary data (e.g. 
the infinite well enclosure). This issue has received some coverage in the literature, both physics-oriented 0-0 and 
purely mathematical 0-0]. See also Ref. [lj for additional references and a discussion of earlier attempts to find 
the spectral solution for the infinite fractional well. 
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In the present paper, the term ”ultrarelativistic” directly stems form the notion of the quasi-relativistic operator 
v 7 —A + to 2 (natural units being presumed) and its mass m —> 0 limit jAj 1 / 2 , see e.g. [lj and (l6| . 


II. THE INFINITE WELL ENCLOSURE: FROM |A | 1/2 TO |A|)/ 2 . 

The Hamiltonian-type expression H = — | A| 1 / 2 + V, with V(x) = 0 for x £ D = (—1,1) C R, is an encoding of the 
Cauchy operator with the Dirichlet boundary conditions (so-called zero exterior condition on R\D) imposed on L 2 (R) 
functions f(x) in the domain of H: f(x) = 0 for \x\ > 1. We point out that the Cauchy operator | A| 1//2 if restricted to 
a domain comprising solely L 2 (R) functions with a support in D and vanishing on R\D is not a self-adjoint operator 
in L 2 (R). 

However, if we consider the action of |A| 1//2 on test functions / £ C^°(D) (infinitely differentiable functions that 

are compactly supported in R ), then the restriction |A|^ 2 of lAp/ 2 to D is interpreted as the Cauchy operator 
with the zero (Dirichlet) exterior condition on R\D and is known to extend to a self-adjoint operator in L 2 (D), fl3| . 
The passage from Cq° (R) to Cq° (D ) ultimately amounts to disregarding any R\D contribution implicit in the formal 
definition (ID- 

Let us discuss the D versus R\D interplay in more detail, by considering the action of | A| 1 / 2 on these C^°(R) 
functions which are actually supported in D, i.e. ip £ C^°(D), while departing from the original nonlocal definition: 

|A|'/y W = -l / t^±V)_m dy , (2) 

n Jr y 



ip{x), x £ (-1,1) 
0, otherwise, 


(3) 


Given x £ (—1,1), we realize that ip(x + y) does not vanish identically if x + y £ (—1,1) i.e. for — 1 — x < y < 1 — x. 
Therefore, the integration ([2]) can be simplified by decomposing R into (—oo < y < — 1 — x) U (—1 — x < y < 
1 — a:) U (1 — x < y < oo). We have: 


|A \ 1/2 iP 


1 

7r 




f 1 x ip(x + y)~ ip( x ) 
J-i-x y 2 


_ 2 ip(x) __ 1 f' 1 x ip(x + y) - ipjx) 

irl-x 2 ir J_ 1 _ x y 2 


(4) 


where the second integral should be understood as the Cauchy principal value with respect to 0, i.e. f\* x = 

linw) [n_ s + /!r]. 

Given x £ (—1,1), let us make a substitution x + y = t in dH) , presuming that now the Cauchy principal value needs 
to be evaluated relative to x. We obtain (note the principal value (p.v.) symbol, introduced in the self-explanatory 
notation) 


,a|./ V = HAA + I I 

7T 1 — 7T J _ 

2ip(x) J x 


1 ip{x) - ip(t) dt = 2 ip{x 


= — lim 

7r e->0 


-1 ( t-x ) 2 
ip{t)dt rl 

x-\-e 


7T 1 — X 2 


-(p.v.) 


ip(x) 


(t - xf 


r - -h«) / 

J x+e {t-x)i * J-1 


t — X 

ip(t)dt 


r 1 ip(t)dt 


(5) 


where [R) refers to the Hadamard regularization of hypersingular integrals (Hadamard finite part, extensively em¬ 
ployed in the engineering literature, [^-[12]). We point out that the troublesome term ^ has been cancelled 
away by its negative coming from the evaluation of (p.v. )[...] in the above. 

__ -1 I c\ 

The second line of the formula (0 can be interpreted as a definition of |A |^3 . The pertinent operator, instead of 
referring merely to C^(D) functions, can be literally applied (extended) to functions ip £ L 2 (D). 

In the literature on the usage of the Hadamard finite part evaluation of hypersingular integrals, it is often mentioned 
that if the (R) integral and the (p.v.) integrals in question do exist, we can relate them as follows: 


|A| D 2 ip(x) 



I —(pv) f 1 ^ dt 
7T dx v J_ x t — x 


( 6 ) 
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We shall employ another version of the Hadamard - Cauchy integral relation, by following a direct integration by 
parts procedure and continually keeping in mind that the involved integrals are (hyper)singular, c.f. also [191 \‘2'A . 

Namely, by invoking the second line of Eq. © and performing integrations by parts before the limit e —> 0 is 
ultimately taken, we end up with: 


|A| D 2 ip{x) 




ip'(t)dt 
t — x ’ 


(7) 


where ip'(t) = dip(t)/dt. 

1 /2 

We are interested in solving an eigenvalue problem |A|^ tp = E ip for the infinite Cauchy well, while interpreted 
in terms of the hypersingular integral equation. For explicit computations we shall employ the Hadamard (finite 
part)-Cauchy (principal value) relation 0: 


E ip(x) H— ( p.v .) 


1,1 


t — x 


= 0. 


( 8 ) 


III. cos(7ra:/2) AND sin(7ra;) ARE NOT EIGENFUNCTIONS OF |A|^ /2 . 


In Ref. 0 we have discussed the validity of counter-arguments against proposed so-far, in the physical literature, 
spectral solutions for Levy-stable infinite well problems [111 [26(. By invoking rigorous mathematical results of EMI 
we have given in na the computer-assisted proofs (elaborated for the infinite Cauchy well) that spectral results 
of [24j-(29| are surely incorrect in the lower part of the spectrum and may be employed at most as approximate 
expressions for higher eigenvalues. In particular a computer-assisted analysis of approximate eigenfunctions shapes 
3 have demonstrated quite clearly that ’’plain” trigonometric functions (like e.g. sine and cosine) are not the 
eigenfunctions for the problem under consideration, see also j26|. 

Presently we shall demonstrate analytically that cos(7ra;/2) is not the ground state function of | A | ^ 2 , so contradicting 
claims of [27H29j . Our method is different from that adopted in Ref [261 ] . The integrations are to be performed in the 
Hadamard sense, c.f. © and that will allow us to introduce basic tools that will be necessary in the subsequent more 
general spectral analysis. 

Let us directly substitute ip(x) = cos(7nr/2) to Eq.®. We shall demonstrate that: 

ahUfl - cihl±£l', 
2 2 J 

(9) 

which surely remains incompatible with any function of the form E cos(7ra;/2), where E > 0 is a constant and 
x £ (—1,1). Here Ci(a;) and Si(a;) are respectively the cosine and sine integral functions, which are defined as follows, 

m : 


11/2 
Id 


nX 1 /ID 

COS — =-(rt) , ,, „ 

2 7r J_ 1 (t — x ) 2 2 


1 cos ^-dt 1 7 tx 

= - cos 


Sihl±£l + sihLA 


1 . nx 
-\— sin — 

2 2 


, f x sin t , 1 sin t , 

Si(a;) = / - dt =-/ - dt 

Jo t ^ J x ^ 


is an entire function on R with the properties Si(oo) = n/2 and Si(— x) = —Si(cc), while 

cos t — 1 


cost , r 

Ci(ai) = — - dt = C + In x + j 

J x t Jo 


t 


-dt 


( 10 ) 


( 11 ) 


is restricted to R + and C = — J 0 °° e _t In t dt = 0.577215665... stands for the Euler-Mascheroni constant. 

We point out that Si (a;) is defined everywhere on R as a continuous and differentiable function. On the contrary, 
Ci(*), x £ R + logarithmically escapes towards —oo as x drops down to 0 Sim. In the vicinity of the well boundaries 
x —> ±1 of D, the logarithmic divergence ln(l — |x|) —> —oo definitely dominates. 

For the direct evaluation of |A|]^ 2 cos with x £ (—1,1) we shall employ the Cauchy principal value formula ©. 
According to Eq. ®, we have: 


I A Id 2 c° s (^-) = 2 ^ p - v -' ) 


1,1 sin(7rf/2) , 1 

V 1 , -dt = -(jp.v.) 


^1—X 


sin[7r(w + a:)/2] 


-l 


t — x 


du = 


' — 1 — X 


(12) 


+ - cos(7nr/2) (p.v.) 


pl—X 

— 1—x 


sm(7ra/2) , 1 . , ... 

- -du + — sui(ttx/2)(p.v.) 

u 2 


r* 1—x 


cos(7tm/2) 


du. 


— 1—x 
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By employing the definition (1101) . we readily get 
rl ~ x sin(7rii/2) 


{p.v.) 


-1— X 


-du = lim 

ej.0 


pi—x - 


-1—x 


sintW^spa+^+gj-rCl-^ 


(13) 


To evaluate (p.v.) f\^ x cos (^ M / 2 ) let us notice that 

f 1 ~ x cos(7tu/2) — 1 I" f~ e f 1 ~ x 1 cos(7tu/2) — 1 

(p.u.) / -cm = lim ' 1 ' 


ej.0 


-1—x •/ e 


du = In + 

1 — x 


/»1+X /»1 — X" 

'o Jo 


cos(7ru/2) 


By employing CD, we get 


, , r 1 x cos(7ru/2) 

(p.u.) / - du = Ci(l — x) — Ci(l + x) 


du. 

(14) 


(15) 


and the identity © readily follows. Clearly, the outcome ([21) is incompatible with E cos(7nr/2), E > 0, as predicted 

in H3, 0-0) 

With all necessary tools in hands, one can easily verify that the the lowest odd would-be (candidate) eigenfunction 
sin(7rx) (that according to [24} ]. 0429} 1 of |A|^ 2 is a faulty guess. Namely, we have: 

|A|^ 2 sin(7ra;) =- (H) j — —ydt = sin(7rx) ^Si[7r(l — a:)] + Si[7r(l + 1 )]^ — cos(7ra;) ^Ci[7r(l — a;)] — Ci[7r(l + z)]^ , 

(16) 

while an expected outcome (according to [13], [2714291 ] ) should be E' sin(7rx), where E' > 0 is a constant. This is 
definitely not the case. We point out that a logarithmic divergence becomes dominant at the boundaries of the interval 
D, that in view of Ci[7r(l — a;)] —> —oo for x f 1 and — Ci[7r(l + a;)] —> +00 as a: J. —1. 

The above discussion easily extends to more general formulas (12311 and (13611 in below, which provide a direct 
analytic demonstration that trigonometric functions of the form cos(n7ra;/2) and sm(mrx/2) with n integer, are not 
the eigenfunctions of |AUj . That invalidates claims to the contrary, appearing in the literature on so-called fractional 
quantum mechanics, [24[, [271429} . 

We note that on formal grounds, the trigonometric functions seem to be valid eigenfunctions if the Fourier multiplier 
representation (c.f. Section I) is ’’blindly” used, ignoring the subtleties related to Fourier integrals of functions with 
support in a bounded domain, SI!- The point is that the primary, mathematically well founded, definition of the 
fractional operator is provided by the integral formula m and not by its Fourier integral version. The latter is merely 
a derived one while on R , 0- If spatial constraints are imposed, we may keep their effects under tight control only 
on the level of Eq. 0, see our considerations of Section II. 


IV. SOLUTION OF THE EIGENVALUE PROBLEM IN THE INFINITE ULTRARELATIVISTIC WELL 

Now we are going to solve the integral equation (]S|). i.e. to deduce the eigenfunctions and eigenvalues of the nonlocal 

operator |A|}/ 2 . We note [§ 3 } that there are no worked out systematic methods (even numerical) of solution of integral 
equations if their kernels are singular, or (that is worse) hypersingular. In below we shall provide an example of a 
successful solution method, based on Fourier series (trigonometric) expansion in L 2 (D). Derivations of approximate 
eigenfunctions and eigenvalues are computer-assisted. The outcomes converge slowly towards ” true” solutions due to 
the singular behavior of Ci at the boundaries of D. 

To find the eigenfunctions and eigenvalues of the nonlocal operator | A|^ , we adopt the following assumptions: 

1. Based on standard quantum mechanical (Laplacian based) infinite well experience and previous attempts, 00 
and 0-0 , to solve the Levy-stable infinite well problem we can safely classify eigenfunctions to be odd or even. 
The oscillation theorem appears here to be valid and the ground state has no nodes (intersections with x axis), 
first excited state has one node, second one has two nodes etc. So, our even states can be labeled by quantum 
numbers k = 0,2,4,6,... while odd states by k =1,3,5,.... 

2. The (Hilbert) state space of the system can be interpreted as a direct sum of odd and even (sub)spaces, equipped 
with basis systems comprising respectively even and odd orthonormal sets of functions in the interval [-1,1]. 



















5 


I/O 

3. In accordance with the infinite well boundary conditions, the function in the domain of (A^ must obey 
i/j(x) = 0 for |x| > 1. In consequence, among various orthogonal sets available in L 2 (D ), we are ultimately left 
with standard trigonometric functions. 

4. The even basis system in L 2 (D) is composed of cosines 


ip k {x) = cos^±i^, J ip k {x)wi{x)dx = 5 kl , k > 0, 


where 5 k i is the Kronecker symbol. For the odd basis system we take the sines 


Xk{x) = sin kirx, J Xk(x)xi( x )dx = Ski, k>l. 


(17) 


(18) 


4. 


We look for eigenfunctions of |A|^ 2 separately in odd and even Hilbert (sub)spaces of L 2 (D). 
Presuming that the Fourier (trigonometric) series converge, for even functions we have 


while for odd functions 


OO 

^e(x) = y a k cos 
k—0 


(2k + l)7rx 
2 


OO 

V>o(x) = E bk sin k'KX. 

k=\ 


(19) 


( 20 ) 


To avoid confusion, we point out that the standard numbering of overall infinite well eigenfunctions begins from 
n — 1 rather then from k = 0 (even case) or k = 1 (odd case) as we have assumed above. We need to have a clear 
discrimination between sine (odd) and cosine (even) Fourier series expansions. The final outcomes will be re-labeled 
according to the traditional lore, i.e. in terms of consecutive integers n = 1,2,.... 


A. Even subspace 

In the present case we substitute the function if> e (x) © into J5]) to obtain 

(2k + l)7rx 


^ a k f k (x) = E ^ a k cos ■ 


/c—0 


/c—0 


where 


f k (x) = -~(U) [ 

n J -i 


1 cos (2fc + 1)7r * 

2 dt = 


1 + 2 k 


sin 


(t-x) 2 
(2k + l)-7rx 


( 21 ) 


( 22 ) 


(2k + l)7r(l — x) ^.(2k+ 1 )tt( 1 + x) 

L 1-—- — V.-1-—- 


+ 


+ cos 


(2k + l)7rx 


Si 


(2k + l)7r(l — x) (2k + l)7r(l + x) 

Dl- 


(23) 


We recall that the functions are singular at x —> ±1 and diverge as ln(l — |x|). Nonetheless, matrix 

elements computed in below prove to be finite. It is the singularity of f k (x) which slows down a convergence of 

1 /2 

approximate expressions for A|^ ip e (x) (finite series expansions of increasing accuracy) to the corresponding ’’true” 
eigenfunctions Eif> e (x). 

Let us multiply both sides of the equation (1211) by +i(x) m and integrate from —1 to 1, while employing the 
orthonormality of ip%(x). The equation (1211) is now replaced by an (infinite) matrix eigenvalue problem 


oo „i 

£ a klki =Ea h + ki = / f k (x)ipi(x)dx , i,k,l = 0,1,2,3,... 

i,fc=0 J ~ 1 


(24) 
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whose solution will be sought for in terms of a sequence of eigenvalue problems for finite n x n matrices, with a 
gradually increasing degree n. 

The set (l24l) is the linear homogeneous system, which, according to Kronecker-Capelli theorem, has a nontrivial 
solution only if its determinant equals zero. This permits to determine the eigenvalues Ek and the coefficients a*, of 
the expansion 03 as the eigenvectors, corresponding to each E^. That separately for each degree n of the involved 
matrix. 

While solving Eq. AMD numerically, the best way to calculate 7 ki is computer-assisted as well (integrals necessary 
to evaluate 7 ^ are no more divergent, so that their numerical calculation is straightforward), but it turns out that 
a number of them can be computed analytically. The analytical calculation permits to establish the fact that the 
matrix (1241) is symmetric, i.e. 7 ^ = 7 ^. and thus the sought for eigenvalues are real. 

In particular, we have 

Ikk = —— + (2k + l)Si[7r(2fc + 1)]. (25) 

7r 

Some exemplary 7 ^ are worth reproducing as well: 


100 = - - + Si W = 1.21531728, 710 = 701 = 6000 - 6 ( 3 ( 3 *) +In 729 = „. 2773259i 

7r 07T 

720 = 702 = (2Ci(7r) - 2Ci(57r) + In 25) = -0.2227035, 

247r 

5 / 15625\ 

721 = 712 = T—- ( 6Ci(37r) - 6Ci(57r) + In " ) = 0.3088509. 

l07r V 729 / 


(26) 


Remark: For the reader’s convenience, let us mention that an analytic evaluation of matrix elements 
can be greatly simplified by taking advantage of worked out indefinite integral formulas (Section 5.3 of 0 ) 
for e.g. f cos(ax)Ci(/3x)dx, f sin(ax)Ci(/3x)dx and analogous integrals with Ci replaced by Si. It is worth¬ 
while to notice that if such integrals contain products of trigonometric functions instead of ’’plain” ones, we 
can always reduce them to one of the listed forms by employing various trigonometric identities. Example: 
2 sin(ax) 0 x 3 ( 72 ;) = sin[(a + 7 ) 2 ;] + sin[(a — 7 ) 2 ;]. 


An explicit form of the matrix (1241) . once we truncate the infinite series at a finite n, reads 



( 700 

710 • 

• InO \ 


710 

711 • 

' 7 nl 

a d = 





\ 7 n 0 

7 nl • 

' Trm / 


(27) 


To find its eigenvalues and eigenvectors we use iterative procedure, considering partial matrices 2 x 2, 3 x 3 etc. The 
eigenvalues of the simplest partial matrix 2 x 2 give the lowest order approximation of ground state and second excited 
state n = 2. The equation for associated eigenvalues reads: 


7oo — E 710 
7 io 7 n - E 


= 0 , 


(28) 


The analytical expressions for Eq and E 2 can be obtained by means of analytical formulas for 7 ^ (|26p . Although com¬ 
putations are cumbersome, one arrives at a reasonable (albeit still far form being sharp) approximation to eigenvalues 
associated with the ground state and first (even) excited state. Using numerical values of 7 ^ (1^51) . we deduce 


E 0 = 1,191256, E 2 = 4.411727 - eigenvalues, (29) 

i/j(E 0 ) = (-0.996257, 0.086437), ip(E 2 ) = (0.086437, 0.996257) - eigenvectors.. (30) 

In other words, the approximate (crude, low order) shapes of the eigenfunctions read 

TTX 37 TX 

i/jq = —0.996257 cos_ ^“ + 0.086437cos —— — ground state, (31) 

7TX 3tTX 

ip 2 = 0.086437 cos — + 0.996257 cos —— second excited state. (32) 

We note here that the reproduced eigenvectors are L 2 (D) normalized, while an overall sign may be negative, which 
is immaterial for the validity of the spectral solution. 
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By increasing the matrix order from 2 to 3, we improve the accuracy with which lowest states are reproduced and 
increase their number by one. We have for eigenenergies 

E 0 = 1.1814891, E 2 = 4.3854565, E 4 = 7.569241. (33) 

It is seen that while one more state appears, numerical outcomes for lowest states are corrected by approximately 1%. 
For the 6x6 matrix we have 

E 0 = 1.1704897, E 2 = 4.35648331, E 4 = 7.52132, E 6 = 10.68291, E s = 13.845025, E w = 17.01393. (34) 

We note that the value E 0 (1341) is quite close to the (still crude) approximate eigenvalue E ga = 37r/8 = 1.1781 deduced 
in Ref. (El E3)- According to [13| the infinite Cauchy well eigenvalues E n become close to (^ — f) — > , as 

n —> oo. Obviously, while passing to higher order matrices the obtained eigen-solutions give better approximations of 
’’true” eigenvalues and eigenvectors in the infinite Cauchy well problem. 

The analysis of numerical values of matrix elements in (1271) shows that these of diagonal elements are much larger 
than the off-diagonal ones. This difference appears to be lowest for joo which equals 1.215, while off-diagonal elements 
take values around 0.3, see (l26l) . For larger k the diagonal elements grow (for example 722 ~ 4.388), while off-diagonal 
values remain close to 0.3. This means that diagonal elements (expression (l25| for even states and (l37l) for odd ones) 
give a fairly good (even if crude) approximation for eigenvalues of the matrix (1271) . Compare e.g. also the first row of 
Table I. 


B. Odd subspace 


We look for eigenfunctions in the form (l20l) . Repeating the same steps as for the even subspace we generate the 
following set of equations 


00 

E b k m=Eb u Vki = g k {x)xi{x)dx, i,k,l = 1,2,3,..., 
i,k =1 J ~ 1 

9k(x) = -- f 

* J -1 


sin knt 

I -1 {t~x ) 2 1 

= fc|sin(/c7ra;) ^Si[fc7r(l — a;)] + Si[/c7r(l + a;)]^ — cos(fc7ra;) ^Ci[/c7r(l — a;)] — Ci[fc7r(l + x)]^ |. 
We find analytically 

Vkk = 2 k Si(2fc7r). 

Eigen-solutions for the 2x2 matrix have the form 

Ei = 2.81019, E 3 = 5.99476 — eigenvalues, 
ip(Ei) = (-0.995891, 0.0905574), ip(E 3 ) = (0.0905574, 0.995891) - eigenvectors.. 


(35) 


(36) 


(37) 


(38) 

(39) 


Two lowest eigenvalues of the 6x6 matrix read Ei = 2.78021, E 3 = 5.93979. In Table I we reproduce the remaining 
four eigenvalues in the 6 x 6 case, in a comparative vein. Namely, we display the computation outcomes for lowest 
six eigenvalues, while gradually increasing the matrix size, from 6 x 6 ,12 x 12, 5000 x 5000 to 10000 x 10000. We 
reintroduce the traditional labeling in terms of i = 1,2,3,4, 5, so that no explicit distinction is made between even 
and odd eigenfunctions. Our results are directly compared with the corresponding data obtained by other methods 
in Refs. (l3L H3 | and (2,0]. 

In Table II we report the change of the ground state energy while increasing the matrix size from 30 x 30 to 
10000 x 10000. It is seen that the third significant digit stabilizes already for 300 x 300 and 400 x 400 matrices. 


C. Graphical comparison 

First, we plot the first four eigenfunctions in Fig. [T] It is seen that qualitatively the states in the Cauchy well at a 
rough graphical resolution level do resemble those (appear to be close) of the ordinary quantum infinite well (deriving 
form the Laplacian). Anyway, we know perfectly (see e.g. Section II) that ’’plain” trigonometric functions, like e.g. 



i 

1 

2 

3 

4 

5 

6 

Diagonal elem. 

1.21531728 

2.83630315 

4.38766562 

5.96864490 

7.53320446 

9.10820377 

Ei6x6 

1.1704897 

2.780209 

4.356483317 

5.9397942 

7.52131594 

9.099426 

Eil2xl2 

1.1644016 

2.7690111 

4.3388792 

5.919976 

7.4952827 

9.0725254 

Enoi x io 4 

1.157791 

2.754795 

4.3168638 

5.892233 

7.460284 

9.032984 

E i( K) [13] Table 2 

1.1577 

2.7547 

4.3168 

5.8921 

7.4601 

9.0328 

Ei(K kms) [Ml Eq. (11.1) 

1.1577738 

2.7547547 

4.3168010 

5.8921474 

7.4601757 

9.0328526 

E i(Z G)M Table VII 

1.1560 

2.7534 

4.3168 

5.8945 

7.4658 

9.0427 

Ei( Z g) [4] Table III 

1.157776 

2.754769 

4.316837 

5.892214 

7.460282 

* 


TABLE I: Comparative table of 6 lowest eigenvalues E % in the Cauchy infinite potential well. Results for matrices of different 
sizes in our approach are compared with spectral data of Refs. 0, m a and [j. First six diagonal elements of the matrix 
(12711 (expressions (12511 and (1371) respectively) are cited for comparison. Note that the numbering of states follows tradition 
(i = 1, 2, 3,4, 5, 6) and refers to consecutive eigenvalues, with no reference to the parity of respective eigenfunctions. 


n (matrix n x n) 

30 

50 

100 

200 

400 

1000 

2000 

5000 

10000 

Eg.s. = Ei 

1.160505 

1.159428 

1.158608 

1.158193 

1.157984 

1.157858 

1.157816 

1.157791 

1.157791 

E2 

2.760953 

2.758572 

2.756705 

2.755742 

2.755252 

2.754954 

2.754855 

2.754795 

2.754795 

e 3 

4.326418 

4.322736 

4.319842 

4.318343 

4.317578 

4.317114 

4.316958 

4.316864 

4.316864 

Ei 

5.904768 

5.900041 

5.896238 

5.894235 

5.893204 

5.892573 

5.892361 

5.892233 

5.892233 

e 5 

7.476052 

7.470114 

7.465334 

7.462812 

7.461511 

7.460714 

7.460446 

7.460284 

7.460284 

e 6 

9.051406 

9.044604 

9.039015 

9.036021 

9.034462 

9.033504 

9.033180 

9.032984 

9.032984 


TABLE II: The matrix n x n-”size evolution” of six lowest eigenvalues of m as n grows. E g . s . stands for ground state energy. 


1 /2 

cos(7rx/2) or sin(7ra;), are not the eigenfunctions of |A|_q . Quite detailed analysis of the eigenfunctions shape issue 
can be found in Ref. 3, where another method of solution of the Cauchy well problem has been tested. 

Since, in the present paper, we employ trigonometric functions as the orthonormal basis system, for low-sized 
matrices (Ell) we deal with visually distinguishable oscillations. These are gradually smoothened with the growth 
of the matrix size. It is instructive to compare approximate shapes of the ground state function, obtained by the 
diagonalization of different-sized matrices. The left panel of Fig. [2]reports the pertinent shapes in case of 3 x 3, 5 x 5 
and 30 x 30 matrices. We note that the qualitative features of the ground state function approximants are practically 
the same for matrices of sizes exceeding 30 x 30. 



_i_i_i__i_i_i_ 

-1 -0.5 0 0.5 1 


X 

FIG. 1: Four lowest eigenfunctions in the infinite Cauchy well, labeled i = 1,2, 3, 4. Outcome of the 10 4 x 10 4 matrix. 
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FIG. 2: Left panel. Comparison of the shapes of ground state functions obtained by the diagonalization of 3x3 (black dashed 
curve), 5x5 (red dash-dot curve) and 30x30 (blue solid curve) matrices. The shape of ground state functions for matrices more 
then 30x30 are identical to that for 30x30. Right panel shows the approximation of ground state wave function (for 700x700 
matrix, solid curve) by the expression (1401) (dashed curve). As both lines are indistinguishable in the scale of the figure, the 
inset depicts the modulus of the point-wise difference of respective curves 


In Ref. [4j an analytic approximation of the ground state function of |A 


1/2 

D 


has been proposed in the form 


4>i{x) = ipg.s.{x ) = 0.921749-^/(1 - x 2 ) cos ax , a = ■ (40) 

In the right panel of Fig. [2] we compare the ground state function (l40l) with that obtained by the diagonalization of 
700 x 700 matrix (which turns out to be close to that obtained by means of the 30 x 30 matrix, see Fig. ?? below). 
It is seen that both functions are indistinguishable within the scale of the figure. The inset in Fig. [2] depicts the 
modulus of the point-wise difference of these functions. Interestingly, although the approximation is non monotonous 
(the difference oscillates), in a large portion of f the interval —1 < x < 1 the difference does not exceed 0.005. 


D. Eigenvalues of |A|^ 2 . 

If compared with the previous methods of solution 0Q and m , our spectral approach seems to be particularly 
powerful if one is interested in the eigenvalues of In fact, we are able to generate an arbitrary number of 

eigenvalues, with a very high accuracy. In Table III we compare several (first 20 and a couple of larger) lowest 
eigenvalues of lA^ and answer how much actually the ’’rough” approximate formula nn/2 — 7r/8 deviates from 
computed E n s. That is motivated by the upper bound formula, 13, 14] (in our notation and for the Cauchy stability 
index a = 1), whose right-hand-side drops down to 0 with n —> oo: \E n — ^ ^ | < 

It is seen from the Table III that although the asymptotic formula delivers pretty good approximation to the 
desirable eigenvalues, the relative error never (except for n = 11) falls below 10~ 3 % as the label number n grows. 
We have actually traced this statement up to n = 500. Moreover, the relative error, as it is seen from the Table III, 
oscillates around 10 -3 %, which means that beginning with n ~ 8 the expression n7r/2 —7r/8 contributes 5 significant 
digits of the ’’true” asymptotic answer. 

Technical comment: We note here that to diagonalize large matrices (30 x 30 and larger) we use the Fortran 
program, based on the LAPACK package. All integrations involved in the evaluation of 'jki and rjki have been 
performed numerically. 
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n 

E n , 5000a;5000 

nir ir 

2 8 

Relative error, % 

Data from [141 

i 

1.157791 

1.178097 

1.75 

1.157773 

2 

2.754795 

2.748894 

0.21 

2.754754 

3 

4.316864 

4.319690 

0.06 

4.316801 

4 

5.892233 

5.890486 

0.03 

5.892147 

5 

7.460284 

7.461283 

0.013 

7.460175 

6 

9.032984 

9.032079 

0.01 

9.032852 

7 

10.602447 

10.602875 

0.004 

10.602293 

8 

12.174295 

12.173672 

0.0051 

12.174118 

9 

13.744308 

13.744468 

0.0012 

13.744109 

10 

15.315777 

15.315264 

0.0033 

15.315554 

11 

16.886062 

16.886061 

5.9-10" 8 

* 

12 

18.457329 

18.456857 

0.0026 

* 

13 

20.027767 

20.027653 

0.00057 

* 

14 

21.598914 

21.598449 

0.0021 

* 

15 

23.169448 

23.169246 

0.00087 

* 

16 

24.740517 

24.740042 

0.0019 

* 

17 

26.311115 

26.310838 

0.0011 

* 

18 

27.882131 

27.881635 

0.0018 

* 

19 

29.452773 

29.452431 

0.0012 

* 

20 

31.023751 

31.023227 

0.0016 

* 

30 

46.731898 

46.731191 

0.0015 

* 

50 

78.148251 

78.147117 

0.0015 

* 

100 

156.689159 

156.686934 

0.0014 

* 


TABLE III: Several lowest eigenvalues of the 5000 x 5000 matrix (1271) are presented. For comparison the approximate formula 
nn/2 — n/8 is depicted together with the relative error \E n — (nn/2 — n/8)\/E n . Independently obtained spectral data (formula 
1.11 in 1 j) are displayed as well. 


V. CONCLUSIONS 

In the present paper we have elaborated a novel, independent from previous proposals, method of an approximate 
solution of the spectral problem of the infinite Cauchy well. Our method is based on the reduction of the initial 
spectral problem for the operator — | A | 1 / 2 3 to the Fredholm-type integral equation with the hypersingular kernel. This 
equation, in turn, can be solved by means of the (Fourier series) expansion with respect the complete set of orthogonal 
functions on the interval — 1 < x < 1 (trigonometric functions which are eigenfunctions of the Laplacian). 

The adopted (Fourier series) expansion method transforms the integral eigenevalue problem © to the eigenvalue 
problem for an infinite matrix. We solve the approximate eigenvalue problems for finite matrices, of the gradually 
increasing size. With the growth of the matrix size, new higher eigenvalues are generated, while lower eigenvalues 
becomes more and more accurate. We demonstrate that the lowest eigenfunctions can be approximately inferred by 
means of the diagonalization of relatively small matrices, like e.g. 30 x 30. We have noticed that the diagonal elements 
of an approximating (finite) matrix give already good approximations for the eigenvalues, see Table I. To obtain the 
eigenvalues with 6 significant digits the diagonalization of matrices of the size 10000 x 10000 of more is necessary. 

The method appears to be a particularly powerful tool to compute the eigenvalues. It can be generalized to other 
fractional (Levy stable) operators, like e.g. —|A|^/ 2 , /x £ (0,2). 
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